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ABSTRACT 

General  diffraction-strain  equations  are  deveJoped  for  the  optical 
diffraction  grating  strain  gages  now  coming  into  use  in  connection  with 
experimental  studies  of  wave  propagation  in  solids.  Equations  are 
derived  for  the  following  cases:  superimposed  gratings  crossed  at  any 
angle  and  subjected  to  arbitrary  surface  strains;  orthogonal  gratings 
subjected  to  arbitrary  rotations;  and  orthogonal  gratings  aligned  with 
principal  strain  directions  and  subjected  to  arbitrary  rotations  during 
strain. 

The  inherent  strain  resolution  capability  of  an  optical  diffraction 
system  is  discussed  and  the  differences  between  this  technique  for  strain 
measurement  and  that  of  x-ray  stress  analysis  are  described. 
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GENERAL  OPTICAL  DIFFRACTION 
STRAIN  RELATIONS 


I.  INTRODUCTION 


The  new  methods  for  analysis  of  mechanical  strains  by  optical  dif¬ 
fraction  techniques  cannot  bt  aiscussed  without  reference  to  the  con¬ 
ceptually  related  method  of  strain  analysis  by  x-ray  diffraction.  One 
might  be  tempted  to  look  to  this  older  method  for  guidance  in  developing 
suitable  equations  describing  diffraction  by  visible  light.  In  fact, 
there  is  very  little  that  can  be  adopted  from  the  established  x-ray 
methods  beyond  the  application  of  the  familiar  Laue  equations  and  the 
liberal  use  of  the  reciprocal  lattice  concept.  The  need  for  the  specific 
equations  newly  developed  in  the  following  report  arises  from  the  desire 
to  measure  strains  in  dynamic  loading  which  are  accompanied  by  surface 
rotations.  These  rotations  can  take  the  form  of  shear  strains  or  of  actual 
rigid  body  rotation.  The  separation  of  strains  and  rigid  body  rotation 
is  imperative  if  the  measurements  are  to  be  meaningful.  Exploring  the 
extent  to  which  this  can  be  accomplished  is  one  of  the  main  goals  of  this 
report . 

Using  the  x-ray  method  (1)  as  a  point  of  departure,  it  should  be 
recognized  that  the  underlying  principles  of  the  optical  method  are  the 
same.  The  differences,  however,  are  significant:  The  lattice  seen  by 
x-rayo  is  a  3-dimensional  array  of  atoms  and  requires  that  3  Laue 
conditions  be  satisfied  simultaneously,  placing  rather  str..ngent  restric¬ 
tions  on  the  phenomenon.  An  attempt  to  measure  a  strain  in  a  single 
crystal  with  the  use  of  monochromatic  x-rays  would  lead  to  a  total  loss 
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of  the  diffraction  spot.  In  the  manmade  optical  grating  formed  into  the 
surface  of  a  specimen  at  least  one  of  the  three  Laue  conditions  is  relaxed. 
Consequently,  a  complete  image  of  a  reciprocal  lattice  section  always 
remains  visible  when  the  monochromatic  light  of  a  laser  is  used.  The  x-ray 
method  requires  polycrystaline  materials  to  satisfy  the  3  Laue  equations 
with  sufficient  frequency  to  produce  useful  information  and  the  results 
are  more  conveniently  expressed  in  the  language  of  Bragg's  law.  Then 
the  two  methods  have  very  little  in  common.  It  might  be  noted  that  the 
only  case  of  an  x-ray  technique  related  to  the  present  work  is  found  in 
the  back-reflection  method  of  orientation  determination  for  single 
crystals  in  which  the  Bragg  reflections  from  the  faces  of  a  crystal  zone 
share  some  rotation  properties  with  the  diffractions  from  one  or  two 
dimensional  optical  gratings.  Comparing  the  two  techniques  in  their 
applications,  the  following  peculiarities  stand  out.  Both  methods 
measure  strains  in  the  surface  of  a  solid;  the  x-ray  method  can  also 
detect  a  strain  perpendicular  to  the  surface  but,  because  of  high 
absorption,  only  in  a  thin  surface  layer.  The  x-ray  method  has  some 
serious  limitations.  Thus,  (a)  the  sample  must  be  crystalline  so  that 
most  plastics  are  automatically  excluded,  (b)  only  elastic  strains  can 
be  detected,  and  (c)  even  with  flash  x-ray  sources,  high  speed  applications 
in  the  microsecond  range  are  not  feasible  yet.  The  x-ray  method  does 
have  an  inherent  strain  resolution  that  is  much  higher  than  that  of  the 
optical  grating.  It  follows  that  the  optical  grating  is  most  suitable 
for  high  speed  measurement  of  plastic  strains  and  is  limited  only  by 
strains  high  enough  to  impair  the  structure  of  the  grating  to  an  extent  that 
it  ceases  to  function  as  an  ordered  diffraction  source. 
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Previous  workers  (2,  3,  4,  5,  6,  7)  utilizing  the  diffraction  grating 
strain  gage  have  derived  their  diffraction-strain  relations  from  the 
equation 

sin0  -  sini  ■  nAm  (1) 

n 

£ 

where  8  ■  the  diffraction  angle  of  the  n  direction  of 

constructive  interference  measured  from  the 
normal  to  the  grating 

i  -  the  angle  of  incidence  of  the  light  measured  from  the  normal 
n  ■  the  order  number 
A  -  the  wave  length  of  the  radiation 
m  -  the  reciprocal  line  spacing  of  the  grating. 

This  equation  applies  only  to  a  line  grating  where  the  direction  of  the 
incident  illumination  is  required  to  be  perpendicular  to  the  lines  of  the 
grating  but  is  not  necessarily  perpendicular  to  the  grating  plane.  When 
this  is  the  case,  the  diffraction  orders  as  seen  on  a  plane  screen  will  lie 
along  a  straight  line.  The  spacing  of  the  orders  along  such  a  line  is 
easily  determined  from  equation  (1) . 

If  the  direction  of  the  incident  light  is  not  perpendicular  to  the 
grating  lines,  the  locus  of  the  diffraction  spots  on  the  observation  screen 
will  be  curved.  When  this  is  the  case  equation  (1)  is  inapplicable.  It 
then  is  necessary  to  use  a  more  general  diffraction  equation  valid  for  an 
arbitrary  direction  of  the  incident  light. 

The  use  of  equation  (1)  involves  some  basic  assumptions  about  the 
behavior  of  the  grating  during  deformation.  If  diffraction-strain  relations 
derived  from  equation  (1)  are  to  be  absolutely  correct  then  the  diffraction 
spots  must  remain  on  a  straight  line  during  deformation  of  the  grating.  For 
this  reason  those  rigid  rotations  of  the  grating  which  would  result  in  a 
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curved  pattern  must  be  excluded.  A  tilt  of  the  grating  plane  about  any 
line  which  ia  parallel  to  the  linea  of  the  grating  would  not  alter  the 
atraight  line  character  of  the  diffraction  pattern.  However,  all  other 
rotatlona  muat  be  either  zero  or  small  enough  to  produce  a  negligible 
effect  on  the  diffraction  pattern.  When  the  grating  linea  are  not 
parallel  to  a  principal  strain  direction,  large  shear  strains  will  rotate 
the  grating  lines.  This  rotation  could  cause  a  curved  locus  of  dif¬ 
fraction  spots  not  adequately  described  by  the  equation.  To  avoid  this 
behavior  the  lines  of  the  grating  are  usually  ruled  in  anticipated 
directions  of  principal  surface  strain.  Thus  a  diffraction-strain  theory 
based  on  equation  (1)  is  limited  to  a  grating  whose  lines  are  parallel  to 
a  principal  strain  direction  and  whose  rotation  axis  is  that  same  principal 
direction.  When  other  rotations  or  shear  strains  are  present  the  use  of 
equation  (1)  could  lead  to  serious  inaccuracy. 

A  disadvantage  of  the  line  grating  strain  gage  is  that  it  measures 
only  one  strain  component.  A  crossed  grating  provides  sufficient  infor¬ 
mation  to  determine  the  entire  state  of  surface  strain.  Unfortunately, 
the  existing  diffraction-strain  equations  can  be  used  with  crossed 
gratings  in  only  very  limited  cases. 

There  are  several  reasons  for  developing  general  diffraction-strain 
relations.  These  equations  can  be  used  to  determine  the  accuracy  of 
existing  ones.  They  apply  to  crossed  gratings  and  exhibit  the  effects 
of  all  rigid  rotations.  In  some  cases  shear  strain  effects  can  also  be 
incorporated  into  the  equations. 

The  purpose  of  this  report  is  to  present  some  general  diffraction- 
strain  equations  which  apply  to  both  line  and  crossed  gratings.  Equations 
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derived  first  are  for  a  crossed  grating  but  with  no  rotations.  These  are 
followed  by  a  description  of  rigid  rotations.  Then  general  diffraction- 
strain  equations  that  include  rotation  effects  are  obtained.  A  discussion 
of  the  resolution  capabilities  of  the  grating  strain  gage  is  also  provided. 


II.  THE  RESPONSE  OF  A  CROSSED  GRATING  TO  ARBITRARY 
SURFACE  STRAINS  IN  THE  ABSENCE  OF  SURFACE  TILT 

A  crossed  grating  has,  as  the  name  implies,  two  sets  of  parallel  lines. 
The  resulting  grid  can  be  approximated  by  a  two  dimensional  doubly  periodic 
lattice.  The  points  of  Intersection  of  the  lattice  lines  can  be  considered 
as  secondary  sources  of  light  in  the  Huygens'  sense.  A  direction  of  con¬ 
structive  Interference  occurs  when  the  light  from  any  four  adjacent 
sources  differs  in  phase  by  integral  numbers  of  wave  lengths.  The  geo¬ 
metry  of  the  undeformed  lattice  is  determined  by  the  lattice  vectors  aQ 

and  b  as  shown  in  Figure  1.  Vector  c  is  normal  to  the  lattice  plane 
0  o 

and  has  the  physical  dimension  of  length  as  do  aQ  and  The  direction 

of  the  incident  light  is  given  by  vector  iQ.  A  typical  direction  of 
constructive  Interference  is  shown  as  vector  d  .  These  vectors  must 


o 

satisfy  the  equations 

(d  (hk)  -  i  )  •  a  -  hX  ,  h  -  o,  +  1,  +  2 . m  (2) 

o  oo  —  — 

(d  (hk)  -  i  )  *  b  -  kX  ,  k  -  o,  +  1,  +  2,  ....  n  (3) 

ooo  —  — 

|dQ|  -  1  (4) 

and  |iQ|  -  1  (5) 


The  subscript  o  refers  to  the  undeformed  configuration.  These  equations 
are  the  well  known  Laue  equations  of  x-ray  diffraction  theory  see,  for 
example,  Azaroff  (8).  The  directions  of  constructive  interference  are 
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completely  determined  if  the  lattice  vectors,  the  direction  of  incidence 

and  the  wave  length  are  specified. 

Equations  (2)  through  (5)  would  also  apply  to  the  deformed  lattice. 

Hence  to  determine  diffraction-strain  equations  one  must  solve  these 

equations  for  both  the  undeformed  and  the  deformed  lattice  vectors. 

A  solution  to  equations  (2)  through  (5)  is  facilitated  by  introducing 

reciprocal  lattice  vectors  a  *,  b  *  and  c  *.  These  vectors  are  determined 

ooo 

from  the  lattice  vectors  by  the  relations 

b  x  c 
0  0 

a  *  »  —  ,r  -  s  (6) 

o  a  •  (b  x  c  ) 

0  0  o 

c  x  a 
o  o 

b  *  «  —  —  — 

o  a  •  (b  x  c  )  (7) 

ooo 

a  x  b 

o  o _ 

c  *  ■  —  ,/r  —  \  (8) 

o  a  *<b  x  c  ) 

0  0  o 

Reciprocal  vectors  (aQ*.  bQ*)  will  be  shown  to  be  directly  related  to  the 
observed  diffraction  pattern.  The  reciprocal  vectors  satisfy  the  following 
conditions 


a  **a  - 

0  0 

1 

a  *«b  ■ 

0  0 

0 

a  *.c 
0  ( 

b  *.a  ■ 
o  o 

0 

b  *.b  - 
0  0 

1 

b  *.c 
o  < 

c  *.a  - 

0 

c  * .b  * 

A  A 

0 

c  *.c 

a  e 

A  solution  is  obtained  by  finding  aQ,  f3  ,  and 'y^  for  which 

d(hk)  ■  a  a  *  +  B  b  *  +7c  *  . 
oo  oo  too 


i  wn — r - ~~ — Tf. 

*1  ■) 


Making  use  of  equation  (10)  in  equations  (2)  and  (3)  gives 


and 


a  ■  i  -a  +  hX 
o  oo 


8  ■  i  *b  +  kX  . 

o  oo 
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(12) 


If  equations  (10)  and  (4)  are  combined  there  results 


1  ■  a  2 1  a  *  | 2  +  8  2|b*|2+'y2|c*|2  +  2a  8  a  *-b 

o'o1  o'o'  f  o  '  a  '  ono  c 


o  '  0 

Equation  (13)  yields 

■\x  1 


c  *  Y  -  + 

1  o  '  'o 


1  -  a  2|a  *|2  -  8  2 lb  *|2  -  2a  8  a  *-b  * 
o  o  o'o  000  o 


(13) 


(14) 


where  the  plus  sign  is  chosen  for  a  reflection  lattice.  One  can  easily 

show  from  equation  (6) ,  for  c  normal  to  a  and  b  ,  that 
n  o  o  o 


a  * 
o 


a  sin(a  ,b  ) 
'o'  o  o 


(15) 


where  Z-(aQ»  bQ)  is  the  angle  between  the  lattice  vectors.  Likewise,  it 
follows  from  equation  (7)  that 


b  * 

1  o  1 


b  sin(a  ,b  ) 
'o'  o  o 


(16) 


From  Figure  2  one  can  see  that 


pb°  ■  la0lsln(VV 


(17) 


and 


P  0  -  |b  | sin(a  ,b  ) 
a  o  oo 


(18) 


where  P,°  and  P  0  are  perpendicular  distances  between  lattice  lines.  The 

reciprocal  vectors  to  aQ  and  bQ  are  also  shown  in  Figure  2.  Reciprocal 

vectors  a  *  and  b  *  define  a  new  lattice  which  is  reciprocal  to  the  lattice 
oo 

of  aQ  and  bQ.  From  equations  (15)  through  (18)  it  follows  that  the  mag¬ 


nitudes  of  the  reciprocal  vectors  are 
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and  |  bQ* |  ■  _1_ 

P  C 
a 


>  (19) 


Equations  (19)  show  the  reciprocal  character  of  the  two  lattices.  The 
angles  between  the  lattice  vectors  of  each  lattice  are  also  related.  From 
equations  (6)  and  (7)  the  scalar  product  of  the  reciprocal  lattice  vectors 


a  *  and  b  *  is 
o  o 


cos(a  ,b  ) 

a  **b  *  - - — — 

o  o  „  o  _  o 


(20) 


P  K 
a  b 

where  equations  (19)  have  been  used  in  the  result.  However,  from  the 
definition  of  the  scalar  product  and  equations  (19) 

C08(a  *,b  *) 


a  *.b  *  «  a  *  b  *  cos(a  *,b  *)  - 
o  o  'o  o  0  0 


o  '  o 


P  °  PK° 

a  b 


(21) 


where  Z_(aQ*,bo*)  is  the  a- gle  between  the  reciprocal  lattice  vectors. 
Combining  equations  (20)  and  (21)  yields 


cos(a  *,b  *)  -  -  cos(a  ,b  ) .  (22) 

oo  oo 

The  solution  to  equation  (22)  that  is  pertinent  to  the  physical  case  is 

£(I  *,b  *)  -  n  -L( I  ,b  )  (23) 

0  0  0  o 

Equation  (23)  relates  the  angle  between  the  real  lattice  vectors  to  the 
angle  between  the  reciprocal  lattice  vectors.  Together  equations  (19)  and 
(23)  determine  the  geometry  of  the  reciprocal  lattice  in  terms  of  the  real 
lattice.  It  will  be  shown  later  that  the  reciprocal  lattice  is  directly 
related  to  the  observed  diffraction  pattern. 

The  solution  for  d  is  found  by  combining  equations  (11),  (12), 

and  (14)  with  equation  (10)  to  obtain 


% 


(hk)  _ 


(i  .a  +  hX)a  *  +  (i  -b  +  kX)b  *  + 

OO  O  00  0 


1  -  (i  .a  +  hX) 2 1  a  *|2  -  (i  .b  +  kX)2|b  *|2  -  2(i  .a  +  hX)(i  .b  +  kX)X 
00  'o'  00  'o'  00  00  ' 


I  - 

—  —  C  * 

a  *,b  *  o 
o  o  — — 

| c  *| 
o 


(24) 


Equation  (24)  gives  the  various  directions  of  constructive  interference  in 
terms  of  lattice  vectors  and  the  vector  of  incidence.  This  equation  is 
used  to  locate  the  orders  on  an  observation  screen  or  the  film  plane  of 
a  camera.  Consider  such  an  observation  screen  as  shown  in  Figure  3. 

For  convenience  the  screen  is  placed  parallel  to  the  lattice  and  the 
direction  of  incidence  is  chosen  normal  to  the  lattice.  The  distance 
between  the  screen  and  the  lattice  is  D. 

For  normal  incidence 


i  .a 
o  o 


0 


and 


i  .b  =0 
o  o 


(25) 

(26) 


Using  equations  (25)  and  (26)  in  equation  (24)  gives 
j  (hk) 

U  —  u/va 

O  O 

.  C 


where  N  T‘ 


hXa  *  +  kXb  * 
o  _  o 

* 
o 


o  i- 


+/  1  -  h2X2|a  *|2  -  k2X2|b  *|2  -  2(hX) (kX)  a  *.b  *  N 
v  o  o  o  o  o 

(28) 


(27) 


is  a  unit  normal  vector  of  the  lattice.  For  normal  incidence  the  direction 
of  zero  phase  difference  as  determined  from  equation  (27)  is 

d  (oo)  -  N  .  (29) 


This  direction  which  is  usually  called  the  zero  order  of  diffraction 
always  corresponds  to  the  direction  of  specular  reflection  of  the  incident 
light.  The  diffraction  spot  on  the  observation  screen  corresponding  to 
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the  zero  order  is  located  by  vector 
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D  (00)  -  Dd  (00)  -  DN 
0  0  0 


where  |d  ^*0  |  is  to  be  determined. 

Vector  D  *hk)  -  D  ^oo)  -  |d  (hk)|d  ^hk)  -  DN^°o) 

A  A  1  A  1  A 


(30) 


The  diffraction  spot  of  order  (h,k)  is  likewise  located  by  vector 

D  ^hk^  -  |d  ^hk^  |  d  ^hk^  (31) 


(32) 


lies  within  the  observation  plane.  It  follows  that 

(D  (hk)  -  D  (oo))  •  N  -  0  . 

0  0  0 

Using  equation  (32)  in  (33)  and  solving  the  resulting  expression  for 

In  (MO 


yields 


|d  (hk) 

1  o 


i  <hk)i 

o  o 


(33) 


(34) 


With  this  result  equation  (32)  becomes 

n  A  (hk) 

D  (hk)  -  D  (°0)  -  ^  ° v  ,  -  DN 

o  o  d  (hk)>N 


(35) 


When  equation  (27)  is  used  in  equation  (31)  there  results 


D  (hk)  -  D  (00) 
o  o 


D  (hXa  *  +  kXb  *) 

-;ck;  ° 

0  0 


where  d  (hk)  .N  -  cos  (®(hk))  - 

O  O  0 


J  1  -  h2x2 1 a0* I  2  -  k2x2 I bo* I  2  -  2 (hX ) (kX)  ao*.bQ* 


and 


@  (hk) 


is  the  angle  between  d 


(hk) 


o  “  o 

equation  (36)  can  be  rewritten  as 


and  the  normal  N  .  Hence 
o 


(36) 


(37) 


....  ,  .  D  (hXa  *  +  kXb  *) 

5  (hk)  .  5  (00) - °° 

0  0  cos©(hk) 


(38) 
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The  vector  DQ^k)  -  D0^k\  within  the  observation  plane,  locates  order 
(hk)  with  respect  to  the  zero  order.  The  factor  1/cos©^^  equation 
(38)  is  a  distortion  factor  that  results  from  using  a  flat  observation 
screen  rather  than  a  spherical  one.  For  diffraction  spots  near  the  zero 
the  value  of  this  factor  is  nearly  unity.  Using  equations  (19)  and  (20) 
in  equation  (37)  gives 


cos 


©(hk)  -  /  l  -  (— )2  -  (— )2  +  2(hX)  (kX) 


- - — t 

;os(a  ,b  ) 


o  o 


a 


P  °  P  ° 
a  b 


Equation  (39)  indicates  that  cos 


©(hk) 


is  nearly  unity  when 


(39) 


and 


hX 

>  0 
b 

kX 


«  1 


<  1 


(40) 


For  those  orders  (h,k)  that  satisfy  inequalities  (40)  the  distortion  in  the 
pattern  is  negligible.  One  can  think  of  a  "distortion"  circle  about  the 
zero  order  within  which  the  distortion  is  negligible.  Outside  of  this 
circle  the  pattern  is  appreciably  distorted.  The  radius  of  the  circle  is 
dictated  by  the  degree  of  distortion  that  can  be  tolerated.  This  is 
illustrated  in  Figure  4  for  an  orthogonal  lattice. 

Within  the  circle  of  Figure  4  equation  (38)  reduces  to 

h  «  0,  +  1,  +  2,  . . . ,  +  m 


D  (hk)  -  D  (oo)  -  DX(ha  *  +  kb  *) 
0  0  0  0 


0,  +1,  +2,  ...,+n 


(41) 


According  to  equation  (41)  the  locations  -  D^00^  of  the  diffraction  spots 

within  the  distortion  circle  are  the  points  of  the  reciprocal  lattice 
(haQ*  +  kbQ*) ,  magnified  by  the  factor  DX. 
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Assume  that  an  orthogonal  lattice  (aQ,  b  )  is  subjected  to  arbitrary 
normal  and  shear  strains.  It  is  assumed  that  the  strains  are  homogeneous 
across  the  lattice  and  that  the  lattice  plane  is  not  tilted  by  the  de¬ 
formation.  The  undeformed  and  deformed  lattices  are  shown  in  Figure  5.  A 
rigid  translation  has  no  effect  on  the  diffraction  pattern  as  seen  in 
the  focal  plane  of  a  well  corrected  lens.  Hence,  the  rigid  translation 
part  of  the  deformation  can  be  ignored  as  long  as  a  lens  is  used  in  the 
formation  of  the  pattern.  The  reciprocal  lattices  of  the  undeformed  and 
deformed  real  lattices  are  shown  in  Figure  6,  The  original  position  of 
order  (h,k)  on  the  observation  screen  is  from  equation  (38) 


\  .  DA  (ha  *  +  kb  *) 

D  (hk)  -  D  <00> - — °-T-  '  S 

0  0  co.©(hk) 

o 

The  position  of  this  diffraction  spot  after  deformation  is 


(38) 


D(hk)  _  D  (oo)  .  DX  (ha*  +  kb»} 

°  Jhk) 

cos0 


,  (hk) 

where:  cos0 


[  ~2  ,  2(hX)(kX)  cos(a,  b) 

Jl-  (jr>  -  (j-T  +  PaP„ 


(42) 


(43) 


0  a 

a*  and  b*  are  reciprocal  vectors  to  the  deformed  lattice,  and  and 

are  the  perpendicular  distances  between  lattice  lines  after  deformation. 

The  vectors  D  _  p  (°°)  ancj  p^k)  _  ^  (00)  are  measured  from  the 

00  0 

diffraction  pattern.  Angles  ^^k)  an<j0(hk)  are  a^so  measurable  quantities. 
They  are  determined  from  the  equations 

tan  ®(hk)  ■  -1^2. 

o  u 

„  .  Hthk)  |D(hk)  -  n  (°o)  1 

and  tan©  -  I _  0  |  .  (45) 


-  (hk)  _  -  (00) 


(44) 


If  a  lens  is  used  then  D  is  replaced  by  the  focal  length  of  the  lens. 


i 


n  &,x  •*-#  ai&& W!*1!  * 


_ t' - 
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The  extensions  of  the  lattice  vectors  are 


a  -  a. 


o'  1  a 1 


and  e.  ■ 

D 


>-1  -  \K\  Id 


-  1  . 


From  equations  (38)  and  (42) 


i ©(h°)  (D  (h0)  -  D  (00)] 

OO  0 

DhX 


,©(ho)  [D(h0)  -  D  (00)] 


For  an  originally  orthogonal  lattice  (aQ,bo)  the  magnitudes  of  aQ*  and  a*  are 


0  P,°  |- 

b  a 


|  a*  |  ■  — —  ■  — - — 

Pb  |a|sin(a,b) 


Combining  equations  (48) ,  (49) ,  and  (50)  gives 
,  coS®(ho)  |d  (h0>  -  d  (oo)  | 

1  O  1  0  ‘ 


and  a* 


coB^ho)  |DCh0)  -do(00> 


| a  I sin(a,b) 

Rewriting  equation  (53)  gives 


sln(a,b)  cos^h>o) 


D  <00) 
0 


However,  from  Figure  6 


£(a,b)  +  £.(a*,b*)  ■  tt  . 


j 


I  W 
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This  equation  also  follows  from  equation  (23)  that  was  developed  for  a 
non-orthogonal  lattice. 

From  equation  (55)  it  follows  that 

sin(a,b)  -  sin(a*,b*).  (56] 

The  extension  of  lattice  vector  aQ  is  obtained  by  combining  equations 
(46),  (52),  (54)  and  (56)  to  get 
cos©(ho)  |D  (h0>  -  D  (oo>  I 

m  o  1  o _ o  1 _  -  1  .  (57] 

a  (h,o)  —(ho)  —  (oo)  .  .  t-.. 

cos0  D  -  Dq  sin(a*,b*) 


A  similar  calculation  gives  the  extension  for  lattice  vector  b  as 

0 

C08®(°k)|D0<0k)  -Do(°o)| 

'  coS|J°k)  |D<ok)  -  Do(oo)|8ln<I*,b*>  "  1  '  <5i 

The  shear  strain  is  given  by 

rab  ’  f  '  '  <5! 

Combining  equations  (55)  and  (59)  gives 

rab  *  L(a*,b*)  -  y  (6( 

where  (_ (a*, b*)  is  measured  directly  from  the  diffraction  pattern.  Due  to 
equations  (44)  and  (45)  cos©^k^  and  cos^^k^  can  t,e  determined  directly 


,©(hk)  . 


T 7|D  (hk)-  D  (00)|2 
+  0  o 


and  cos/ 


TT |D  (hk)  -  d  (00)|2 
0 


Hence  to  compute  the  extension  of  lattice  vector  aQ  one  need  only  measure 
three  quantities.  These  are: 


v.  nirtC 1  .wW*/1  V  iW  *  <  .4 


*fw  t 


,.**»mjt'a»iW  "  ■  ^-Y*' ' 
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^(a*,b*)  , 


and 


|D  (h6)  -  D  (°0)|. 

-L>  5 ?°°> 

D(h0)  Do  . 


To  compute  the  extension  of  lattice  vector  bQ  two  additional  measurements 


are  needed.  These  are: 


Id  <ok)  -  d  (oo) 


and  |D(ok)  -  D  (oo) 
o 


The  necessary  quantities  to  be  measured  are  shown  in  Figure  7.  Once 
these  quantities  are  measured  the  entire  state  of  strain  can  be  determined. 


I 


III.  THE  DESCRIPTION  OF  RIGID  ROTATIONS  OF  A  CROSSED  GRATING 

The  most  general  diffraction-strain  equations  should  include  the  effects 
of  an  arbitrary  finite  rigid  rotation  of  the  grating.  Since  finite  rota¬ 
tions  cannot  be  treated  as  vector  quantities  it  is  first  necessary  to 
select  coordinates  which  locate  uniquely  the  angular  orientation  of  the 
grating.  These  coordinates  should  provide  the  position  of  the  rotated 
grating  with  respect  to  its  original  position. 

Suppose  that  an  orthogonal  lattice  (ao,bQ)  suffers  an  arbitrary  rigid 
rotation.  Then  the  final  position  of  the  lattice  can  be  specified  by  a 
set  of  Euler  angles  6,0)  as  shown  in  Figure  8.  The  rotated  lattice 

vectors  are  (a.b)  as  shown.  Vectors  c  and  c  are  the  normals  to  the 

o 

lattice  before  and  after  the  rotation.  Angle  6  is  the  angle  of  tilt  of 
the  lattice  plane  about  line  u. 

Angle  0  establishes  the  axis  of  tilt  of  the  grating  which  lies  in 
the  grating  plane.  The  total  angle  of  rotation  of  the  lattice  about  the 
normal  is  0+0.  It  is  emphasized  that  the  Euler  angles  determine  the 


f 


t 


} 


16 

final  orientation  of  the  lattice  without  regard  to  the  actual  aequence  of 
rotations  that  produces  the  final  configuration. 

The  rotation  can  also  be  described  by  an  orthogonal  matrix  A  that 
relates  the  unit  lattice  vectors  of  the  two  configurations  as  follows 


ui  ‘  Vj  ‘'j_-  1i2-3  (63> 

where  e^  are  original  unit  lattice  vectors  and  u^  are  final  ones.  Since 

matrix  A  and  the  Euler  angles  both  describe  the  same  rotation  they  must 

be  related.  To  determine  this  relation  it  is  convenient  to  take  the 

rotation  as  three  separate  rotations.  The  first  of  these  is  a  rotation 

about  the  normal  c  of  the  grating  as  shown  in  Figure  1.  The  unit  lattice 

o 

vectors  of  the  rotated  lattice  are 

ei  "  aijej  i,i  "  1,2,3  (64> 

where 


a  ■ 


cosl^/  ainljj  0 
-sinl jj  cos \fj  0 
0  0  1 


(65) 


The  second  rotation  is  a  tilt  of  the  plane  of  the  grating  about  line  u>  of 
Figure  9  as  shown  in  Figure  10.  This  rotation  produces  unit  lattice 
vectors  e''  as  given  by 

e-  -  b  ij  i*j  -  1,2,3  (67) 

where 


■  M 


10  0 

0  cos0  sin0  .  (68) 

0  -sin0  cos0 

The  final  rotation  is  about  normal  c  of  the  grating  as  shown  in  Figure  11. 


iu  i 


mi  ,o 


% 


I 


<  ,V  -p- 


Unit  vectors  of  this  configuration  are 


ui "  cijej 


i* j  -  1,2,3 


where 


c  ■ 


cob(£)  sinC^  0 
-sin (p  cosCj!)  0 
_  0  0  1 
Combining  equations  (64),  (67)  and  (69)  gives 


N 


ui  ■  cijbjkaktet 
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(69) 


(70) 


i,j,k,t  -  1,2,3  (71) 

In  matrix  notation  equation  (71)  becomes 

u  -  cbae  (72) 

where  u  and  e  are  unit  column  vectors.  If  matrix  A  is  defined  as 

A  ”  cba 

°r  AU  ■  Ctpbp,aqj  l.J.P.q  -  1.2.3  (73) 
then  equation  (63)  follows.  Since  the  rotation  matrix  A  is  the  product  of 
the  three  successive  rotation  matrices,  it  follows  that 

cos<f>co6^-8i«fecos0siny  cos<j)sin^+sin<j)cos0cosyA  sin4>sin0 
■  '-sin<i>co8y*-cos<£cos0sin^.  -sin^>sin^+cos4)cos0cos^  coscj>sin0  .  (74) 

8in0sin^  -sin0cosl/r  cos0 

To  produce  the  final  configuration  of  the  grating  three  success  Ive  rota¬ 
tions  were  used.  However,  the  actual  rotation  of  the  grating  into  this 
final  configuration  is  arbitrary.  Any  combination  of  rotations  that 
produce  the  same  final  configuration  will  yield  the  same  Euler  angles  and 
matrix  A.  In  fact,  the  final  orientation  can  be  reached  by  two  rotations 
in  at  least  two  different  ways.  The  total  rotation  can  be  achieved  by 
first  a  tilt  of  the  grating  plane  followed  by  a  single  rotation  about 


the  grating  normal.  If  the  plane  is  tilted  through  angle  0  about  line  uj 


7XXT-. 


the  unit  lattice  vectors  become 


<  '  aijej 


1,J  -  1,2,3 


where 


“2  2 

cos  y.+cos0sin  p  sin^cosyi(l-cos0)  -sinOsin^. 

2  2 

sin^cos^>(l-cos0)  sin  ^+cos0cos  p  sin0cos^. 


N  ■ 


1 


i  - 
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(75) 


(76) 


sin@sinjt»  -sinOcos^  cos0 

This  is  followed  by  a  rotation  (</>+$<,)  about  the  new  grating  normal.  The 
last  rotation  produces  the  unit  lattice  vectors 

ui  ■  Vi  1,j  ■ 1>2>3 

where 


(77) 


cos  ((4+0)  sin(^+0) 

i 

-sin(^+0)  cos(v*f0)  0  j  ,  (78) 

0  0  1  j 

p  is  the  angle  that  locates  line  u>,  and <p  is  measured  from  line  u>  as  before. 
The  combined  matrix  of  these  two  rotations  is  given  by 

A  -  b'a'  (79) 

where  A  is  given  by  equation  (74),  a'  by  equation  (76)  and  b'  by  equation 
(78). 

The  final  configuration  can  also  be  reached  by  first  a  rotation  about 
the  normal  of  angle  (<fi+pO  followed  by  a  rotation  about  line  oj  of  angle  0. 
The  first  rotation  is  characterized  by  matrix 

cos(^+<p)  sin(M*0)  0 

-sin(^,+0)  cos(^+<p)  0  .  (80) 

0  0  i 

This  one  is  followed  by  a  rotation  about  line  o>  with  matrix 


t  *.  \ 


i 


cos  0+cos0sin  0  sin0cos^ (cos0-l) 

2  2 

sin0cos/(cos0-l)  sin  <p+ cos0cos  cp 
-sin&sinf,  -sinBcos/- 

The  total  rotation  matrix  is  then 


sin0sint/. 
sin0cos0  | . 
cosfi 
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(81) 


A  =  c'b'  (82) 

where  A  is  also  given  by  equation  (74). 

Regardless  of  the  number  of  rotations  used  to  achieve  the  final 
grating  orientation  the  values  of  0,  0,  and  0  are  the  same.  These  angles 
are  measures  of  the  net  rotation  resulting  from  any  rotation  history. 

Angle  0  is  the  resultant  angle  of  tilt  of  the  grating  plane.  It  is  the 
angle  between  the  normal  before  rotation  and  the  normal  after  rotation. 

Angle  (0+</>)  is  the  sum  rotation  of  the  lattice  about  the  grating  normal. 

This  angle  is  measured  about  the  moving  normal  as  it  is  carried  through 
the  rotation  process.  Angle  0,  establishes  the  axis  of  the  resultant  tilt 
of  the  grating  and  is  designated  as  line  u>  in  Figures  8,  9,  10,  and  11. 

The  choice  of  Euler's  angles  to  describe  grating  rotations  facilitates 
the  analysis  of  the  diffraction  problem.  It  will  be  shown  that  these  angles 
can  be  measured  directly  from  the  diffraction  pattern.  In  fact,  angles 
0  and  Qcan  be  measured  by  observing  only  the  zero  order  of  diffraction. 

To  determine  angle  0  it  is  necessary  to  record  the  movement  of  some  of 
the  diffraction  spots  surrounding  the  zero  order. 


IV.  THE  RESPONSE  OF  AN  ORTHOGONAL  GRATING 
ALIGNED  WITH  PRINCIPAL  STRAIN  DIRECTIONS 
AND  SUBJECTED  TO  ARBITRARY  RIGID  ROTATIONS 
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o 


Diffraction-strain  equations  will  be  derived  for  an  orthogonal  lattice 
aligned  with  principal  strain  directions  and  subjected  to  arbitrary  surface 
rotations.  For  moderate  angles  of  tilt  of  the  grating  it  is  shown  that 
the  effect  of  a  surface  tilt  is  almost  exactly  equivalent  to  that  of  an 
appropriate  change  in  the  direction  of  the  incident  light.  This  is  not  the 
case  for  large  angles  of  tilt  since  the  plane  of  observation  is  fixed 
relative  to  the  incident  beam. 

Consider  an  orthogonal  lattice  (a^b^)  as  shown  in  Figure  12. 

Vector  c  is  normal  to  the  lattice.  The  unit  lattice  vectors  are  e, ,  e„, 
o  12 

and  e^.  The  incident  illumination  is  normal  to  the  grating  before 
deformation.  Directions  of  constructive  interference  are  determined  from 
the  equations: 


d  (ht)  .1  -  hA 

O  0 


d  (hk).b  =  kA 
o  o 


h  =  0 ,  +1 ,  +2 ,  . . .  +  m 
k  =  0 ,  +1 ,  +2 ,  . . .  +  n 


where 


d 


(hk) 


(83) 

(84) 

(85) 


The  subscript  o  refers  to  the  original  unstrained  lattice.  Equations  (83) 
and  (84)  can  be  rewritten  in  terms  of  the  unit  lattice  vectors  as 


'  (hk)  "  hA 

a  .e.  =  - 

o  1  i  — 

a_ 


and 


j  (hk)  ; 
do  ,e2 


kA 


(83) ’ 

(84) ' 


A  solution  for  d  ^k^  is  found  by  determining  a .  ^kk^  that  satisfy  the 
o  1 


relation 


d  (hk)  = 


(hk)" 

al  el  * 


(hk); 


f  (hk); 

a3  e3 


(86) 
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It  is  assumed  that  the  lines  of  the  crossed  grating  are  ruled  in  the 
anticipated  directions  of  the  principal  surface  strains.  When  this  is  the 
case  the  deformation  will  preserve  the  orthogonality  of  the  grating.  In 
addition  to  the  normal  extensions  the  lattice  suffers  arbitrary  rigid 
rotations.  The  directions  of  constructive  interference  at  any  time  during 
the  deformation  are  determined  by  the  equations 


(d(hk) 

-  V 

hX 

•  u  =  - 

h  ■  0,  +1,  . 

•  •  W| 

(92) 

1  3  1 

(d(hk) 

-  V 

kX 

2  Ib| 

k  ■  0,  +1,  . 

. .  +  n, 

(93) 

ld(hk) 

1  - 1  , 

(94) 

ig  - 

1  , 

(95) 

22 


where  iQ  is  a  unit  vector  in  the  fixed  direction  of  the  incident  light, 
and  U2  are  the  current  unit  lattice  vectors,  and  |a|  and  |b|  are  the 
magnitudes  of  the  respective  lattice  vectors.  Figure  13  shows  the  deformed 

*  A  A  *  /U1_\ 

lattice  referred  to  the  original  unit  lattice  (e^,  e2,  e3).  Vector  dv  ' 
defines  the  direction  of  order  (hk)  as  the  grating  is  being  deformed.  To 
solve  equations  (92)  through  (95)  it  is  assumed  that 


d°'k)  .  B^  -  B^  +  B2(hk)u2  +  B3(hk)u3 
(hk) 


(96) 


where  B^  are  to  be  determined.  Using  equation  (96)  in  equations  (92) 
through  (94)  there  results 


.  (hk)  ;  A  ,  hX 

B,  -  i  *u  +  —  , 
1  01  11 

1*1 

(hk)  ;  A  .  kX 

J2  "  VU2  +  —  * 


and  |do^kk^|2  ■  1  ■  [ B, (hk) ] 2  +  [ B2 (hk) ] 2  +  [ B3 (hk) ] 2 
Solving  equation  (99)  for  B3^k^  yields 


83 (hk)  -  +J  i  -  [B1(hk)]2  -  [B2(hk)]2  . 

Combining  equations  (96)  through  (100)  there  results 

a(hk)  -  (i  .;1  +  — )“i +  d  -2  +  -^)a2  + 

0  1  HI  1  0  2  |b|  2 


J 


.  f4  A  1  hX.  2  '  ,  kX .  2 ' 

1  -  (i  *u.  + - )  -  (i  -u,  +  - )  u,  . 

0  1  in  02  ibi  3 


(97) 

(98) 

(99) 


(100) 


(101) 


Equation  (101)  provides  the  directions  of  constructive  interference  during 
deformation  in  terms  of  the  rotated  unit  lattice  vectors. 

The  plane  of  observation  is  placed  parallel  to  the  original  grating 
plane  at  a  distance  D  away.  When  a  lens  is  used  in  the  formation  of  the 
diffraction  pattern,  the  observation  plane  corresponds  to  the  focal  plane 
of  the  lens.  In  this  case  D  is  replaced  by  the  focal  length.  Vector  D 


s-.ifsw# «» ?  ■ 
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to  the  observaclon  plane  locaCes  Che  diffracCion  spoc  of  order  (hk) . 


This  veccor  Is  relaced  Co  d 


(hk) 


by 


n:  (hk) 

n  <hk>  .  ^2 
0 


Dd  (hk) 
o 


,  (hk)  - 
do  "e3 


COS 


0(hk)  * 


(102) 


where 


0(hk) 


Is  Che  angle  becween  d 


(hk) 


and  e^.  The  observaclon  screen  and 


vecCors  D^00^  and  DQ^kk^  are  shown  in  Figure  3.  From  equaClon  (38)  ic 


follows  ChaC 


D  *hk)  -  D  (oo)  -  D  e  +  e  ) 

0  0  in  1  ib  i  2 


cos 


0(hk) 


(103) 


where 


co,©(hk)  -  Il  -  ^i2 


“o|2 


2  2 
k  A 

l*ol2 


(104) 


This  vecCor  lies  in  Che  observaclon  plane  and  joins  orders  (hk)  and  (oo) . 

As  a  resulc  of  deformacion  Che  diffracCion  spocs  will  move  Co  new  locaCions 
defined  by  vecCor  D^kk^ 


-(hk)  ^  Dd 


This  veccor  is  relaCed  Co  d^kk^  by 
(hk) 


;<hk)  : 
a  .  e. 


(105) 


where  d^kk^  is  given  in  equaCion  (101).  IC  is  necessary  Co  find  d^kk^  in 
Cerms  of  Che  original  unic  laCCice  vecCors  e^,  e2»  e^.  The  unic  laCCice 
vecCors  before  and  during  deformacion  are  relaCed  by 


-  A^e^  i» j  -  1,2,3  (106) 

where  A  is  Che  roCaCion  maCrix  given  in  equaCion  (74).  When  equaCion  (106) 
is  used  in  equaCion  (96)  Chere  resulcs 

*(hk)  „  Bi(hk)Ui  -  Bi(hk)Aijej  .  i,j  -  1,2,3  (107) 
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This  can  be  rewritten  as 


:(hk)  ..(hk); 

d  ei 


(108) 


where 


fl.(hk)  (hk). 

6j  Bi  Au  • 


(109) 


If  equation  (108)  is  used  in  equation  (105)  the  result  is 


D<hk) 


De;(hk)ej 


,  - (hk) 


(110) 


Now  vector  D^kk^  -  D  will  lie  within  the  observation  plane  as  shown  in 

o 

Figure  14.  Due  to  equation  (110)  this  vector  is 

(hk);  .  .(hk);  ft.(hk);  . 

D(hk)  _  -  (oo)  „  ^ _ fi  _  ;  .  _ ;i_+„g2 _ V  (111) 

O  ..(hk)  De3  ..(hk)  •  uil) 

e3  33 


Equation  (111)  can  be  written  as 


n  (°°)\  au- 

C0B®  (D _ Do  )  -  (g'(hk)ei  +  S2(hk)e2) 


.(hk) 


(112) 


where  8-j  ■  cos0  and  ©'  is  the  angle  between  order  (hk)  during 

deformation  and  the  original  zero  order  before  deformation.  From  equation 
(112)  the  new  position  of  the  zero  order  is 


(°°) /n(°°)  n  (oo). 
COS0  (D  -  Dq  ) 


(e;(00,e1  +  62(00)e2)  . 


(113) 


Equation  (103)  is  rewritten  as 


,©<!>«  d  (hk)  -  d  (oo)  -  (Ja- :  Jsi.  . 

°  -2 - 5  - -  Ia0l  1  |bj  2 


The  factor  cos 


©(hk) 


(114) 


in  equation  (114)  removes  the  distortion  in  vector 


(D  ^kk^  -  D  ^00^)  which  is  caused  by  a  plane  observation  screen.  Sub- 
o  o 

tracting  equation  (113)  from  (112)  gives 


l 
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c,.,,ihk)  <D(hk)  *  D  (00))  -  cos  (00)  (D<00)  -  D  (oo)) 

0  _ 0 _  ■ 

D  D 

[BJ(hk)  -  B£(oo)]e1  +  [02(hk)  -  ^2 (°o) ] e2  (115) 

The  cosine  factors  In  equation  (115)  can  also  be  thought  of  as  distortion 
factors.  From  equation  (109) 


(hk)  (hk) 

61  6i  Ail 

i  -  1,2,3 

(116) 

and 

(hk)  (hk) 

B2  Pi  Ai2 

i  -  1,2,3 

(117) 

If  equations  (116)  and  (117)  are  used  in  equation  (115)  the  result  Is 

cos Abk)  (D(hk)  -  Do(oo))  _  cosioo)  (D(oo)  -  D  (oo>) 

D  6 


(6 


(hk)  -  6i<oo>)Ail;1  +  (6i(hk) 


-  81(00))Ai2e2  . 


1  -  1,2,3 


(118) 


Due  to  equations  (97),  (98)  and  (100)  it  follows  that 


6, (hk>  -  6 


(oo)  hX 


a 


(hk)  (oo)  kX 
p9  -  p9  ■  - 

2  2  ibi 


(119) 

(120) 


and 


(hk)  (oo) 
S3  63 


[ \  ,  hX  ,2  ~  ,  kX  x  2 

/I  -  (i_ >u.  + - )  -  (1  .u,  + - )  - 

J  O  1  |  —|  O  2  |T*| 


/ 


<5  A  A  A  /»  ft 

1  -  (1  .U  )  -  (1  .uj 

O  1  0  l 


(121) 


Expanding  equation  (118)  and  making  use  of  equations  (119)  through  (121) 
in  the  result  gives 

COS0  (D  -  Do  )  COS0  (D  -  DQ  ) 

5  D  " 

^  (cos4>cos^-  siivf>co80sin{[,)  +•  (-sinc^cos^-cos^cosQsin^) 

|  M  |b| 

- '  / - 


1  r*  “  .  hX  .2  “  “  .  kX  .2  j  -  *2  “  “  .2 

1  -  (i  .u.  +  - )  -  (i  .u,  +  - )  -  /  1  -  (i  -u.)  -  (i  .u0) 

0  1  ii  0  2  ibi  J  01  02 


X 


r 


sin^sin© 


el  + 
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(cos^sin^+sity^cosQcosjfc)  +  (-sin£sin^+  cos#>cos0cos&) 

Jal  lb! 


■  I  .  hA  .2  *  *  .  kA  ,2  /,  “  ,2  *  *  .2  .  „ 

+l- ll  ~  (io‘Ul  +  ^")  "  (io*U2  +  J  1  ~  (W  ‘  (i"‘U^  I* 


V“2J 


(-sin0cos^)J  e 2  • 

Since  the  light  is  originally  normal  to  the  grating 


(122) 


i 

o 


It  follows  that 

VUi  “  e3,Ui  "  ®3  ,Aij ej  "  ~Ai3 * 
From  equation  (124)  and  equation  (74) 


i  -  1,2,3 


(123) 


(124) 


i^u^  m-A^3  -  -  sin^sinB 

A  A 

and  i  ,u„  «  -A--  -  -  cosd>sin0 

o  i  C  J 

Using  equations  (125)  and  (126)  in  equation  (121)  gives 


(125) 

(126) 


(hk)  (oo) 
3  "  63 


l  -  (-^-  -  sin0sin8)2  -  (-^-  -  cos^sinS)2  - 

J  HI  |b| 


j  2  jp 

/ 1  -  (sin0sin8)  -  (cos^sinS)  (127) 

For  sufficiently  small  values  of  sine  the  value  of  g^°°)  will 

be  small  provided  that  the  following  inequalities  are  satisfied  for  order 

(hk) : 


JhA 

HI 

kA 


<<  1 


«  1 


(128) 


These  inequalities  are  satisfied  provided  that  the  line  spacing  is  large 
compared  to  A  and  if  those  orders  near  the  zero  are  used. 

Since  Jl  -  x  *  1  -  \  (129) 

for  x  small  compared  to  1,  equation  (127)  is  approximated  by 


:|/v  *«?» 


6  (Ilk)  _  6  (oo)  .j^ine—  +  cos^sine  —  -  h  X  ,  -  k  -■■ , 

3  3  HI  HI  2|I|2  2|b|2 

Using  equation  (130)  in  equation  (122)  yields 
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(130) 


co^hk)  (D<hk>  -  Do(00)) 

- 


COS0  (D  -  Dq  ) 

- 


(  cos0cos^-  sin0cos0sin^)+  (-sir^cos^;-  cos0cos0sin^) 

M  |b| 

-  2  2  2  2 

,  ,  ,  ,  „hA  .  .  ,  „kA  hr  k  X  s  , 

+  (sin0sin0 -  +  cosAsin© - - - - »)  sinj6sin0 

|I|  |b|  2|a|2  2|b|2 

— -  (cos#>sin^,+  sin0cos0cos^)  +  (-sinpsin^f  cos<pco80co8tt) 

r»i  ibi 

2  2  2  2 

+  (simfsine— —  +  coBdrainS— - - —  ■' — r  -  k  -  ,)  (-sinBcoafc)  |  e,  (131) 

HI  HI  2HI 2  2  H| 2 

Equation  (131)  simplifies  when  the  angle  of  tilt  0  is  small  enough  so  that 

sin0  #  0 
and  COS0  r*  1. 

Equation  (132)  are  satisfied  with  less  than  one  percent  error  for 
values  of  0  less  than  eight  degrees  (0.14  radians). 

For  this  case  equation  (131)  reduces  to 

(hk)  /„(hk)  -  (oo),  „„oJ°0)  /n(°0)  n  (°°)\  - 

CO80  (D  -  Dq  )  -  COS0  (D  -  Dq  )  - 

5  5 


(132) 


cos(0+^)  -  sin($+^,)  e.  + 

-3-  sin(£-ty)  +  cos^+ty.) 

|a|  1  b  |  J 

1 a  1  1  b  | 

L  o  .2,2 

L  .2,2  J 

e2  (133) 


where  terms  that  contain  0 


HI2 


and 


Ibi2 


are  omitted. 


Equation  (133)  can  be  used  to  obtain  diffraction  strain  relations  for 
angles  of  tilt  that  satisfy  equations  (132)  and  orders  (hk)  that  satisfy 
inequalities  (128) .  These  assumptions  would  lead  to  negligible  error  in 
computing  the  position  of  the  diffraction  spot  of  the  second  order  from 
a  grating  of  5000  lines  per  inch  with  a  surface  tilt  of  eight  degrees. 


'  #«F 


***** 
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In  using  equation  (133)  one  is  not  restricted  to  small  angles  of  rotation 
about  the  grating  normal.  However,  the  angle  of  tilt  must  satisfy  equation 
(132). 


The  extensions  of  lattice  vectors  a  and  b  are 

o  o 


-  1*1  -  !*0i  'll  -  i 


•V 

Ibl  -  lb 


and  e.  ■ 

D 


W-- 1. 


i*j 


(134) 


(135) 


The  values  of  | |  and  |bQ|  are  determined  from  the  original  diffraction 
pattern.  From  equation  (114) 


1  cos 


0(ho)  ip  (ho)  _  p  (oo) 


and 


i  DhX 

o' 

i  co.©(ok)  |n  (ok)  -  n  <00)  | 

— —  ■  o  o  o 


where  cos©^ho^  - 
o 


DkX 

1 


fT  (ho)  _  p  (oo)  I  2 

O  0 


(136) 


(137) 


(138) 


and  cos 


@(ok) 


/ 


Id  (ok)  -  d  (oo)|2 
1+0  0 


(139) 


The  value  of  |a|  is  determined  from  equation  (133)  by  considering  order  (ho). 
From  equation  (133) 


cos 


© 


(ho)  ,=r(ho)  „  (oo) 


(Dv,,w/  -  D 


l  «>•$  '  0> 


(oo)  ,^(oo)  _  p  (oo) 


hX_ 

a 


el  +  hX  sin{£+^)  e2 


(140) 


.  A  V*V*.*0*-J 
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The  individual  terms  on  the  left  hand  side  of  equation  (140)  are  determined 
from  measurements  of  the  spacing  of  the  diffraction  pattern.  Taking  the 
magnitude  of  both  sides  of  equation  (140)  and  solving  the  resulting  expres¬ 
sion  for  | a |  gives 

j  a  |  ■ _ DhA _ .  (141) 

l 

cos@  <D  "  D0  >  "  cos@  (n  -  D0  > 


A  similar  calculation  gives 

Ibl  * 


DkA 


•  (142) 


(ok)  /-(ok)  —  (oo). 

o 


COS0  (D  -  Do  )  -  COS0  (D  -  D 


(oo)  ,-(oo)  -  (oo) 

o 


Using  equations  (136),  (137),  (141),  and  (142)  in  equations  (134)  and 
(135)  yields  the  lattice  extensions 


cos 


@(ho)  i-  (ho)  _  -  (oo) 


(h°)  /n(h°)  n  (°°)\  r„J°°)  ,n(°o)  n  (°°)\ 

COS©  "  D0  )  “  COS0  (D  -  D  ) 


-  1  (143) 


and  e,  = 

D 


COS 


0(ok)  |-  (ok)  _  -  (oo) 


-  1 


(ok)  ,-(ok)  -  (oo). 


(oo)  /Tr(oo)  -  (oo) 


cOS(m)  (D  "  Dq  )  “  COS(m)  (D  -  D 


(144) 


The  distortion  factors  appearing  in  the  denominators  of  equations  (143) 
and  (144)  are  determined  from  measurements  of  the  spacing  of  diffraction 
spots  together  with  the  equation 


(hk) 

cos© 


|D<hk)  -  D  (00)|2 
1  +  o 


(145) 


D 

Equation  (122)  would  predict  the  position  of  the  diffraction  spots 
for  large  angles  of  tilt  of  the  grating  plane.  For  angles  of  tilt  that 


■ 


1 
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satisfy  equations  (132)  very  little  error  results  in  using  the  simplified 
equation  (133) . 

The  Euler  angles  of  the  surface  rotation  are  determined  by  measurements 


of  the  spacing  of  the  diffraction  spots  together  with  the  equations 

(D(00)  -  D0(00))1 
“■^■-(5(00)  _5°«.o)  ’ 


(146) 


tan  26  ■ 


|-(oo)  _  -  (oo) 

' _ 0 

D 


(147) 


,  .  .  ...  (ho)  (ho)  —  (oo).  (oo)  (oo)  r-  (oo) 

and  tan  ty+tf>)  -  |_CO80  (D  -  PQV  ')  -  cos^  (nv  -  PQ  (!48) 

r„D  Jh°)  /7:(ho)  -  (oo).  (oo)  ,~(oo)  n  (°°)J 

Icos0  (D  -  Dq  )  -  COS0  (D  -  Dq  )J  x 

where  the  subscripts  1  and  2  denote  the  respective  components  of  the  vector 
bearing  the  subscript.  Equations  (146)  and  (148)  together  with  Figure  15 
indicate  that  angles  ^  and  0  can  be  measured  by  observing  the  movement  of 
only  the  zero  diffraction  order.  To  determine  angle  <f>  it  is  necessary  to 
observe  at  least  one  other  diffraction  spot. 

It  should  be  noted  that  one  could  construct  a  net  having  the  general 
appearance  of  Figure  4  that  could  be  used  in  conjunction  with  a  recorded 
diffraction  pattern  to  measure  the  rotation  angles  of  the  grating.  This 
finely  divided  net  would  be  drawn  for  a  fixed  grating  to  screen  distance 
or  equivalently  a  fixed  focal  length.  The  physical  dimensions  of  the  net 
coordinates  would  be  degrees  of  arc  or  radians.  Such  a  net  is  shown  in 
Figure  16  and  corresponds  to  the  Greninger  chart  for  x-ray  diffraction. 

In  using  the  net  one  would  simply  lay  a  transparency  of  the  recorded 
diffraction  pattern  over  the  net.  By  aligning  a  row  of  zero  diffraction 
orders  (h,o)  or  (o,k)  with  the  horizontal  curves  of  the  net  one  could 


,J.y 
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determine  the  orientation  of  each  lattice  vector  a  and  b,  following  the 
method  used  to  find  the  orientation  of  a  zone  axis  in  x-ray  crystallography 
(1) .  The  angles  involved  in  this  method  are  not  identical  with  the  Euler 
angles  used  in  the  foregoing  discussion. 

It  can  be  shown  that  equation  (133),  which  was  obtained  from  equation 
(131)  by  neglecting  higher  order  terms,  is  the  exact  solution  to  a  different 
diffraction  problem.  This  equivalent  problem  corresponds  to  having  a  zero 
tilt  of  the  grating  plane  but  a  change  in  the  direction  of  incidence  so 
as  to  produce  the  same  zero  diffraction  spot  on  the  observation  screen 
as  would  result  from  a  tilt  of  the  grating.  The  equivalent  problems  are 
shown  in  Figure  17.  In  this  figure  the  plane  of  incidence  is  shown  as  the 
plane  of  the  paper.  The  positions  of  the  various  diffraction  spots  on  the 
observation  screen  for  the  equivalent  problems  are  very  nearly  the  same  as 
long  as  the  angles  of  tilt  are  small.  For  large  angles  of  tilt  the  two 
problems  are  not  equivalent.  The  diffraction  spots  would  not  occupy  the 
same  positions  when  the  tilt  is  large.  In  this  case  equation  (133) 
cannot  be  used  since  the  higher  order  terms  are  no  longer  negligible.  It 
should  be  noted  that  the  reason  the  two  problems  are  not  equivalent  for 
large  angles  of  tilt  is  the  use  of  a  fixed  observation  screen  or  camera 
position.  If  the  observation  screen  were  rotated  in  such  a  fashion  as  to 
remain  parallel  to  the  grating  as  it  is  tilted  then  the  two  problems  would 
be  identical.  However,  this  would  be  impossible  to  do  during  a  dynamic 
experiment  of  short  duration. 


32 


V.  STRAIN  RESOLUTION  OF  A  DIFFRACTION  GRATING 


The  Rayleigh  criterion  will  be  used  to  find  a  theoretical  estimate 
of  the  strain  resolving  power  of  an  optical  diffraction  grating.  According 
to  the  Rayleigh  criterion  two  diffraction  spots  are  barely  resolved  when  the 
intensity  maximum  of  one  falls  on  the  first  zero  intensity  point  of  the 
other.  The  situation  is  illustrated  in  Figure  18.  In  the  figure  A0  is  the 
half  angular  width  of  the  diffraction  order.  This  angular  width  is  given 
by 


A0 


Xm 


Ncos0 


(149) 


where  X  *  wave  length 

m  *  reciprocal  line  spacing 

N  -  total  number  of  lines  in  the  grating 

til 

0^  ■  diffraction  angle  of  the  n  order. 

For  a  derivation  of  equation  (149)  see  Longhurst  (9) .  This  equation  is  not 
restricted  to  normal  incidence.  The  shift  in  the  n^  diffraction  order 
caused  by  a  strain  £  is 


A6 


-nXme 

COS0 


(150) 


where  n  is  the  order  number.  For  a  discussion  of  equation  (150)  see 
Douglas  et  al .  (6).  The  smallest  shift  resolvable  is  the  value  given  by 
equation  (149) .  Hence  it  follows  that 

Xm  nXme 


Ncos0  cos0 
n  n 


(151) 


By  solving  equation  (151)  for  | e |  one  obtains  an  expression  for  the 
smallest  strain  that  can  be  resolved.  This  result  is 


I  e  I  ■ 

min 


_1 

Nn 


(152) 
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The  strain  resolving  power  is  then  defined  as 

R  -  Nn.  (153) 

e 

Equations  (152)  and  (153)  predict  the  inherent  ability  of  the  nth  dif¬ 
fraction  order  to  resolve  a  strain.  The  chromatic  resolving  power  for 
wave  length  X  and  order  n  is 

X / AX  ■  R^  ■  Nn.  (154) 

Thus  the  grating  can  resolve  a  strain  or  change  in  wave  length  with 
equal  ability.  The  total  number  of  lines  in  the  diffraction  grating 
strain  gage  is 


N  -  L  m 
g 

(155) 

where  L 

g 

is  the  length  of  the  gage. 

Using  this  result  equations  (152) 

and  (153) 

become 

Id  -  1 

1  1  min  nmL 

g 

(156) 

and 

R  -  nmL  . 

e  g 

(157) 

For  example,  a  grating  with  reciprocal  spacing  10,000  lines/in.  and  length 
of  0.1  in.  can  resolve  a  strain  of  0.1%  in  the  first  order.  The  corresponding 
resolving  power  of  this  grating  would  be  1000  for  the  first  order  of  dif¬ 
fraction. 

It  is  emphasized  that  the  resolution  as  obtained  by  the  Rayleigh  crit¬ 
erion  is  a  theorteical  estimate  only.  Due  to  variations  in  line  spacing 
real  gratings  will  not  actually  resolve  the  strains  predicted  by  equation 
(156) .  Grating  irregularities  caused  by  large  strains  also  contribute  to 
a  decrease  in  resolving  power.  Hence,  resolution  values  obtained  by  the 
Rayleish  criterion  should  be  taken  as  limiting  values  only. 
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VI.  CONCLUSIONS 

Diffraction-strain  equations  were  first  developed  for  a  crossed 
grating  subjected  to  arbitrary  surface  strains  in  the  absence  of  surface 
tilt.  From  the  derived  equations  it  is  seen  that  one  need  only  observe 
the  displacements  of  two  active  orders  to  determine  the  entire  state  of 
surface  strain,  since  the  zero  order  would  not  be  displaced  in  the  absence 
of  surface  tilt.  A  direct  measurement  of  the  shear  angle  can  be  made  from 
a  record  of  the  diffraction  pattern. 

Rigid  rotations  of  the  grating  plane  were  considered.  Fuler  angles 
were  introduced  as  coordinates  of  the  angular  orientation  of  the  grating. 
These  angles  were  shown  to  be  Independent  of  the  actual  sequence  of 
rotations  involved,  but  were  directly  related  to  the  tilt  of  the  grating 
and  the  rotation  about  the  grating  normal.  Equations  were  provided  for 
determining  the  Euler  angles  from  measurements  of  the  spacing  of  the  dif¬ 
fraction  pattern. 

Diffraction-strain  equations  were  derived  for  an  orthogonal  grating 
aligned  with  principal  strain  directions  including  the  effects  of  rigid 
rotations.  The  effect  of  a  surface  tilt  was  found  to  be  the  same  as  a 
change  in  the  direction  of  incidence  for  small  angles  of  tilt.  For 
large  angles  of  tilt  and  a  fixed  observation  screen  this  would  not  be 


the  case. 
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